#include <math.h>

//矩孔初始化
void init_rec(double *temp , int Im_W , int Im_H , int hol_w)
{
	int i,j;
	for(i=0 ; i<Im_H ; i++){
		for(j=0 ; j<Im_W ; j++) *(temp+i*Im_W*2+j*2+0)=*(temp+i*Im_W*2+j*2+1)=0.0;
	}
	for(i=Im_H/2-hol_w/2+1 ; i<Im_H/2+hol_w/2 ; i++){
		for(j=Im_W/2-hol_w/2+1 ; j<Im_W/2+hol_w/2 ; j++) *(temp+i*Im_W*2+j*2+0)=255.0;
	}
}
//圆孔初始化
void init_circ(double *temp , int Im_W , int Im_H , int R)
{
	int i,j;
	for(i=0 ; i<Im_H ; i++){
		for(j=0 ; j<Im_W ; j++) *(temp+i*Im_W*2+j*2+0)=*(temp+i*Im_W*2+j*2+1)=0.0;
	}
	for(i=Im_H/2-R+1 ; i<Im_H/2+R ; i++){
		for(j=Im_H/2-R+1 ; j<Im_H/2+R ; j++){
			if((sqrt((i-Im_H/2.0)*(i-Im_H/2.0)+(j-Im_W/2.0)*(j-Im_W/2.0))) <= R-1)
			*(temp+i*Im_W*2+j*2+0)=255.0;
		}
	}
}